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Abstract 


Radiation  and  scattering  problems  are  formulated  as  optimal 
control  problems  in  which  either  a  current  or  surface  impedance 
is  sought  from  a  class  of  admissable  functions  which  optimizes  a 
functional  of  the  scattered  far  field.  In  both  cases  the  existence 
of  an  optimal  solution  is  proven.  In  the  linear  (radiation)  case 
constructive  algorithms  for  finding  the  optimal  solution  are 
presented . 
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1.  luTRQDUCTIOH 

In  this  paper,  we  review  results  obtained  by  the  author 
in  close  collaboration  with  R.Z.  Kleinman  [  l],  i_  2],  |”  3j,  [4]. 
The  class  of  problems  referred  to  in  the  title  arise  in  the 
study  of  time  harmonic  electormagnetic  and  acoustic  radiation 
and  scattering.  As  such,  these  problems  are  problems  in  exterior 
domains  in  contrast  to  the  more  common  problems  with  bounded 
domains  usually  encountered  in  the  literature  of  distributed 
parameter  control  systems.  Such  problems  require,  in  addition 
to  the  usual  boundary  condition,  a  so-called  "radiation 
condition"  in  order  to  guarantee  that  the  problem  has  a  unique 
solution.  This  condition  at  infinity  can  then  be  used,  as  in 
classical  potential  theory,  to  recast  the  problem  which  is 
originally  stated  in  differential  form,  into  a  problem  Involving 
integral  operators  acting  on  appropriate  spaces  of  functions 
defined  on  the  boundary  of  the  scattering  obsticle  or  radiating 
structure.  It  is  this  formulation  of  the  control  problems  in 
terms  of  boundary  integral  equations  that  we  will  use  in  our 
discussion  of  the  scattering  problem. 

First,  however,  we  will  consider  a  simpler  linear  radiation 
problem  with  Dirlchlet  conditions  under  the  assumption  that  we 
have  available  a  Green’s  function  appropriate  to  the  exterior 


domain.  This  will  help  to  focus  the  nature  of  the  boundary 
control  problems  we  consider  and  will  allow  us  to  mention 
some  numerical  results  for  the  control  of  radiation  patterns. 

In  the  last  section,  we  will  turn  to  the  more  complicated 
non  linear  optimization  problem  with  a  Robin  boundary  condition 
and  without  the  assumption  that  the  appropriate  Green's  function 
is  available.  In  our  treatment  of  this  non  linear  problem 
of  optimal  control,  we  will  use  a  pair  of  boundary  integral 
equations  derived  in  u  3  ].  the  unique  solution  to  which  affords 
a  solution  of  the  exterior  boundary  value  problem.  That  we 
deal  with  a  pair  of  boundary  integral  equations  rather  than 

a  single  equation  is  a  result  of  the  well-known  non- unique ness 

/  '  ^ 

of  solutions  of  the  first  of  these  equations  at  interior 
eigenvalues  of  the  homogeneous  Dirichlet  problem.  It  must 
be  emphasized  that  these  exceptional  values  of  the  wave  number 
do  not  correspond  to  exceptional  values  for  the  exterlon 
scattering  problem  but  are  the  result  of  the  particular 
reformulation  of  the  exterior  boundary  value  problem  in  terms 
of  Integral  equations.  Various  means  have  been  proposed  to 
circumvent  this  difficulty  (see  e.g.  [6]  and  [7]).  Here, 
we  follow  [ 8  ]  by  Introducing  a  second  equation  involving  the 
normal  derivative  of  the  double  layer  potential  which  restores 
uniqueness  for  real  wave  numbers. 

2.  fCP.HULATIOH  C£  TH-  LIniAR  PR03LLP 

Let  P  be  a  simple  closed  ourve,  sufficiently  smooth  (CT 

2 

is  more  than  sufficient  for  our  purposes)  in  i.  ,  let  A  be  the 
unbounded  region  determined  by  T  and  let  A*  denote  the 


complementary  region.  In  particular,  the  curve  P  will  have  a 
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a  unit  normal  n(q.)  which  varies  cont  inuously  on  P  .  ie 
emphasize  that  we  shall  take  n(q)  to  be  the  outward  drawn  normal 
with  respect  to  Sit  .  Further,  we  whall  write  V*  n"  and  *>/^n* 
to  denote  the  normal  derivative  when  P -4 ptP  from  ilc  and  SI 
respectively. 

■ie  choose  the  origin  of  a  rectangular  coordinate  system 
in  the  interior  of  Si*  and  let  (r,  9)  be  polar  coordinates  of 
and  arbitrary  point  Pi  n .  The  three  classical  boundary  value 
problems,  Dirlchlet  (D) ,  Neumann  (il),  and  3obln  ( 3.)  problems 
for  the  Helmholtz  equation  in  the  exterior  domain  SL  consist  of 
finding  a  function  u  ,  defined  in  SI  ,  satisfying  (1,1)  and  (l,iii) 
below  as  well  as  the  appropriate  boundary  condition  in  (1,11)  : 


(1) 

(V‘ 

k2)  u  *  0 

in  SL  ; 

f  (D) 

u  »  h 

(il) . 

*  u/  *  n  *  h 

on  r  ; 

1(H) 

lu/ >  n  u  ■ 

h 

(ill) 

lim 

J  \  >  u/  *  r  - 

.2 

iku|  do  *  0 

3r 


where  is  a  circle  of  radius  r  lying  entirely  in  jfl  .  This 
last  condition  is  the  radiation  condition  referred  to  above. 

For  the  problem  of  Dirichlet  (D)  the  Green's  function 
has  the  form 

(2)  j(?,Q)  -  (1/4)  H(*U|?-Ql)  +  g(?,4) 

where  g(?,«;/  (  C  (SI)  and  is  the  Hankel  function  of  the 

first  kind,  the  fundamental  solution  of  the  Helmholtz  equation 
consistent  with  the  radiation  condition.  Since  P  is  bounded 
one  may  employ  the  asymptotic  properties  of  G  for  large  r 
viz.  G ( P , <; )  ■  (e1,kr/r*)  $(  •  ,-t)  +o(l/r)  to  represent  the 
solution  u  in  the  far  field  as 
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(3)  u  (  j?  ;  h )  =  (elkr/r^ 


r^)  j  h 


(q)  Van  >;(6  ,q)  ir  +  o(i/r). 

q.  a 


The  regularity  of  G  allows  us  to  define  the  far  field  pattern 
2 

?l  L  '0,2ir  )  in  terms  of  the  Hilbert-ochmidt  integral  operator 
K  :  L2(  P  )— *  L2(0,2<r  )  generated  by  the  kernel  d  n  C-(  •  ,q) 


(4) 


f(  0) 


- 

Jr 


h(q)  ^ ^  n^  G(  8  ,q)  dPq. 


Two  of  the  problems  commonly  discussed  in  the  engineering 

literature  of  antenna  theory  fit  naturally  into  the  present 

context.  de  refer  to  them  as  the  identlf ic?.tlon  problem 

and  the  synthesis  problem.  In  the  former,  we  are  given  an 

actual  far  field  pattern  and  are  asked  to  find  the  boundary 

data  h  whose  associated  solution  has  f  as  its  far  field. 

In  the  latter,  we  are  given  a  desired  far  field  pattern  g 

and  are  asked  to  find  the  boundary  data  h  producing  a  far 

2 

field  closest  to  the  desired  pattern  in  the  L  sense.  In  teres 
of  the  integral  operator  k  Introduced  above,  the  first  of 
these  is  the  problem  of  solving  the  integral  equation  of  the 
first  kind  kh  =  f  with  f  given  in  the  range  of  H,  while 
the  second  is  that  of  finding  an  absolute  minimum  for  the  form 
l\  .di  -  g  i)  .  3oth  of  these  problems  fall  into  the  class  of 
problems  called  ill-posed  in  the  sense  of  Hadamard.  The 
Interested  reader  may  refer  to  the  paper  of  Angell  and  dashed  [5 
for  a  complete  discussion  of  these  problems  as  they  arise  in 
antenna  theory. 

Here,  it  is  our  purpose  to  present  a  different  optimization 


problem  in  which  the  quantity  to  be  optimized  is  the  power 
radiated  into  a  specified  portion  of  the  far  field.  Specifically, 
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Specifically,  we  define  the  far  field  power  in  a  measurable 
set  «f  C  i. 0,2  if]  by 

(2tr  2  r  2  ? 

(5)  ^(h)  :=  *(9)  lf(0)|  d*  =  *(0)  \Khp(0)  d  9 

*  J  o  Jo 

where  d(^)  is  the  characteristic  function  of  the  set  . 

Having  introduced  this  quadratic  functional,  we  nay  pose  a 

meaningful  optimization  problem. 

2 

Let  UCL  (P)  be  a  closed,  bounded,  convex  set. 

We  will  refer  to  the  elements  of  U  as  admissible  controls. 

The  optimal  control  problem  is  now  to  find  an  hQ£  U  giving 
the  functional  ^(h)  its  maximum  value  i.o.  find  hQt  V 
such  that  =<  (hQ)  \  Q^Ch)  ,  for  all  h  t  U. 

Since  the  form  d,(hj  =  (  kh,kh) ,  p  is  both 

*  l  (o,2«r) 

convex  and  weakly  continuous,  it  is  easy  to  prove  the  following. 

2 

Theorem  JL:  If  UCL  (0,2  tr)  is  closed,  bounded  and  convex 
then  there  exists  an  element  hQt  bdy  U  which  is  optimal 
with  respect  to  the  performance  index  Q*  defined  by  (5). 

We  remark  that  this  statement  not  only  asserts  the 
existence  of  an  optimal  solution,  but  also  gives  an  additional 
condition  which  leads  to  a  "bang-bang"  result  for  certain 
special  choices  of  the  class  U  of  admissible  controls 
(see  the  appendix  in  1.1  ]  for  details).  The  requirements  th't 
the  set  J  be  a  closed,  convex  and  bounded  set  in  L  (P  )  still 
allows  enough  flexibility  to  allow  us  to  model  a  number  of 
constraints  suggested  by  physical  or  design  considerations. 

In  the  engineering  literature,  the  former  are  usually  referred 
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to  as  realizability  conditions ,  and  Include  constraint  sets  of 
the  following  forms: 

(a)  U  =  (h|  L^IP  )  \  |h_U$  1  \  (inputs  with  bounded  energy); 

(b)  u  =  \h  t  a1(r  )  1  Uh«<l  }  (  inputs  with  bounded  energy 

and  oscillation) ; 

(c)  U  =  |  h  «  Lc(  P  )  \  ^(q)ih(qj(A^(q),  a.e.  ,  ^ ^  c  C(P  ; 
while  design  considerations  nay  lead  to  control  domains  defined. 

by  state  constraints .  e .3. 

(d)  U  =  )|  |(  1  -  ot  (0)  )  f (0)1  $  --J  (bounded  siae  lobes 

he  remark  that  it  is  case  (c)  for  which  we  can  show  that  there 
exists  a  maximizing  sequence  of  "bang-bang"  controls  since  the 
extreme  points  of  the  set  U  in  (c)  are  functions  of  the  form 

Vo  Y  +Vi^  1  tihere  is  a  measurable  partition 

of  P  .  7or  full  proofs  the  reader  is  referred  to  1]. 

2 

In  the  case  that  :J  is  the  unit  sphere  in  L  (p),  the 
observation  that  the  functional  can  be  expressed  in  terms 

of  the  symmetric  operator  R:=  k-'ar  a,  where  X*  is  the  adjoint 
of  a,  leaas  to  the  following  result. 

Theorem  2:  If  % ^  is  the  largest  eigenvalue  of  the  operator 

?.  and  h  is  a  corresoonding  normalize-’  e i.renvector ,  then 
0  ~  ” 

sup  -i.(h)  =  sup  Uh,h)  2(n)  »  Xa  ■  *iAh  j. 

*n|«  1  fchu=  1  h  (  1  0  *  0 

Cne  may  now  Implement  constructive  methods  to  compute 

tne  optimal  control  by  considering  a  family  of  subspaces 
2 

of  L  (p),  ultimately  cense  in  this  function  space,  and  con¬ 
sidering  tne  finite  dimensional  problems  ?  (Xu  -  Xu)  =  0, 

n 

n  =  1,2,...,  w.nere  P  is  the  projection  associated  with  the 


subspace  A  .  Standard  estimates  for  the  rate  of  co:ive rgence 
n 

are  available  e.g.  one  has  the  estimate 
(«>  -  pn,3“ 

where  X  is  the  eigenvalue  for  the  finite  dimensional 

o,n 

problem  and  A  is  the  true  eigenvalue  (see  Hikhlin  i  9  "!)  . 

o  ■  “ 

This  method  has  bee  implemented  on  the  Burrough's  6R00  at  the 
University  of  Delaware  using  standard  13 UL  eigenvector-eigenvalue 
routines  to  produce  highly  directed  single  and  dougle  beam 
patterns  for  a  circular  antenna  of  radius  a  for  various 
values  of  ka  from  1  to  20. 

!•  IMS.  -f ONLINE AR  PP-03LE1! 

In  this  section,  we  will  consider,  not  a  radiation 

problem,  but  a  problem  in  scattering  wherein  we  assume  that  an 

2 

object,  again  with  a  sufficiently  regular  boundary,  in  E  , 
Interacts  with  a  known  incident  field,  u.1  ,  to  produce  a 
scattered  field  us.  In  this  case,  us  is  a  solution  of  the 
Helmholtz  equation  in  the  unbounded  region  SI  which  satisfies 
the  radiation  condition,  while  u*  is  assumed  to  be  a  solution 
of  the  Helmholtz  equation  in  J.1  . 

The  nonlinear  character  of  the  problem  arises  from  the 
form  of  the  boundary  condition.  Under  suitable  restrictions 
on  the  geometry  of  the  scatterer  and  constitutive  parameters, 
among  which  is  the  requirement  that  the  radius  of  curvature  be 
large  relative  to  the  skin  depth  (see  e.g.  [_io])  the  transition 
conditions  at  the  surface  of  an  imperfectly  conducting  scatterer 
may  be  replaced  by  what  is  usually  called  an  impedence  boundary 
condition.  Here,  such  a  cond.ltion  takes  the  form  of  the 
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0* 

Tobin  boundary  condition  Su/'bn  +  9  u  =  0,  where  nL  (P) 
represents  the  surface  impedence.  It  is  this  surface  Impedence 
which  we  ta.:e  as  the  control  and,  by  requesting  that  e*  lie  in 
a  preass igned  closed,  bounded  and  convex  (and  hence  weak*- 
sequentially  compact)  subset  of  L  (r),  we  x:\  -  again  consider 
the  constrained  optimization  problem  of  finding  that  admissible 
control  whicn  affords  an  absolute  maximum  to  the  functional 
k  m  defined  by  ( 5) . 

In  this  problem,  the  expression  for  the  far  field  depends 
explicitly  on  the  product  of  the  control  and  the  scattered  field 
and  it  is  here  that  the  nonlinearity  arises.  Cur  treatment  is 
more  complicated  than  that  followed  in  the  previous  case  not  only 
due  to  this  nonlinearity  but  also  in  that  we  do  not  assume  that 
we  i'ave  explicit  knowledge  of  the  Green’s  function  for  the 
region.  Instead,  we  rely  on  the  Helmholtz  representation , 
knowing  the  fundamental  solution  If ( ?,« )  =(-i/2)  H ^  (k  |  P-i|  )  , 

to  derive  boundary  integral  equations  whose  unique  solution  may 
then  be  used  to  construct  the  scattered  field  in  the  exterior 
domain  £1  . 

The  scattered  and  Incident  fields  being  regular  solutions 
of  the  Helmholtz  equation  In  .£1  and  respectively,  we  may 
write 

(7)  us  =  3 (  du'Van)  -  D( us )  and  u1  =  DU1)  -  3(luVi  n) 
where  3  and  D  are  the  usual  single  and  double  layer  operators 
with  kernels  >f(P,q)  and  (P,q)  respectively.  If  K 

denotes  the  Integral  operator  with  kernel  n  (p,CJ  ,  then 

is  the  integral  operator  for  D  and  the  usual  jump  relations 
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for  the  normal  derivative  of  the  single  layer  and  for  the 
double  layer  with  density^*,  namely 

)  (p)  =  (I-:-  X)^i(p) 

■*An'(o»i!'p)  =  (-1  +  X)  u.(p) 

(8)  p  r  _  r 

lim  (Du)  (P)  =  (-1  +  X*)u  (p) 

P-+P  1  ' 

^li m_(Dyi)(P)  =  (I  +  I>y*.  (p) 

together  with  the  relation  (7)  and  the  boundary  condition 
lead  easily  to  a  pair  of  integral  equations  which  may  be 
written  as 

(I  +  3r  +  K*)u  =  2u^ 

(9)  _  . 

(~r I  +  K*<r  +  D  )u  =  2iuVdn 

n 

where  the  operator  D  is  defined  by 

(10)  (D nfi)(p)  :=  V&flp  f  j/iil)  ^^^(p.q)  . 

T 

Jot.e  that  the  integral  operators  all  depend  on  the  wave  number 
1:  although  we  have  not  explicitly  indicated  such  dependence. 

The  values  of  the  wave  number  for  which  there  exist  non-trivial 
solutions  of  w  +  Xw  =  0  will  be  called  characteristic  values 
of  ( -a)  and  a  classical  result  is  that  k  is  such  a  character¬ 
istic  value  if  and  only  if  it  is  an  eigenvalue  of  the  interior 
homogeneous  Dirichlet  problem.  As  described  above,  while  the 
exterior  Robin  problem  has  a  unique  solution,  at  interior 
eigenvalues  of  the  Dirichlet  problem  the  first  of  these  boundary 
integral  equations  has  a  multiplicity  of  solutions  and  it  is  for 
this  reason  that  we  consider  the  pair  of  equations  (9). 

The  analysis  in  !.  3  ]  chows  that,  for  appropriate  choices 
of  k,  and  in  particular  for  k  real,  the  system  (9)  has  a 


10 


unique  solution  which  can  then  be  used,  to  construct  a  solution 
of  the  original  scattering  problem.  Specifically,  we  have 
the  following  result  !.  3  J: 

Theorem  Let  »  £  L*“  (D  and  suppose  that  k  satisfies  the 
conditions  Im  k  fcO,  Ia(k»)  *>0.  Then  there  exists  a  unique 
solution  u£L2(D  of  the  pair  of  boundary  integral  equations 
(9).  Moreover,  the  function  u  =  u  +  u  is  a  solution  of  the 
exterior  Robin  problem 

(i)  usc2<n),  u,  -»u/an  t  L2(P )  ; 

(ii)  (72  +  k2)  us  =  0  in  SL  , 

2  2  i 

(V  +  k  )  u  =0  in  «£X  ; 

(iii)  lira  [  |l>us/jr  -  ikuJ  |  d3  =  0  ; 

r-**  J3r  r 

(iv)  "bu/8n  u  a  0  a.  s.  on  P  j 
if  and  only  if 


uS  =  ?3(-'iui/an  -  cr u)  -  fD(u  -  u1) 
where  u  is  the  unique  solution  of  the  system  of  boundary 
integral  equations  (9). 


Returning  to  the  optimization  problem,  it  is  possible 
to  represent  the  far  scattered  field  in  the  form 
(11)  f(  •  )  -  K  (  ^  uVfcn  +o-u)  -  K2us 

2  2 

where  and  X2  are  compact  operators  from  L  (P)  to  L  (0,2*W). 
Using  the  cost  functional  given  in  (5).  we  can  pose  the 
following  optimization  problem: 

Given  UCL* (f )  a  closed,  bounded,  convex  set,  find 
«r  0 1  U  for  which  (  O'  Q )  %  i.  ^  ( or  )  for  all  «r  rU. 

Mote  that  the  results  described  above  guarantee  that, 
given  et'J  there  exists  a  unique  solution  \iia)  of  the  system  (9). 


-  11  - 

If  we  let  A  be  the  set  of  all  pairs  (a1,  u(«*)}  ,  rf  U,  we 
may  prove  the  following  theorems  (see  f_  4  ]  for  complete 
details) . 

«•  2 

Theorem  4:  The  set  AC  1  (p) *  L  (P)  is  bounded  in  the  product 
topology. 

Theorem  The  set  A  is  closed  relative  to  we  ale* -converge  nee 
of  the  «r  and  strong  convergence  of  the  u  . 

The  ore  m  6:  The  map  <r  *— *f  is  continuous  from  the  weak*-topology 
of  L  (P)  to  the  strong  topology  of  L'  (0,2  vr). 

Theorem  £:  There  exists  an  element  <r  t  U  (and  consequently 
a  pair  ^Q,u(  <r  Q )  j  f  A  )  such  that  Q^(  <r  Q)  2  ( cr  )  for  alia*  c 


Clearly,  this  last  theorem  follows  immediately  from 
Theorem  6  if  one  tabes  into  account  the  form  of  the  cost 
functional.  The  techniques  used  to  establish  Theorems  4,  5, 
and  6  are  similar  and  we  conclude  our  discussion  with  a  sketch 
of  the  proof  of  Theorem  4  in  the  hope  that  it  will  give  the 
reader  sufficient  insight. 

Proof  of  Theorem  4:  Suppose  that  the  set  A  were  not 

bounded.  Then,  since  U  is  bounded,  we  may  choose  a  sequence 

\  (a  .  u  )1  C  A  for  which  the  set  {cr  I  is  bounded  but  \lu  1I-* 
mm*  *  m  m 

Let  N*  :=  u  /  a  u  ||  ,  and  choose  a  subsequence  for  which 
1  m  m  "  ra" 

W*  |  W  i  W  2 

(relabeling)  qr  -  ■»  tr  ,  T  »  and  ^  m  ■  ■  — »  p  L  (P) 

Then  the  functions  ''f  satisfy  the  integral  equation 
(12)  (I  +  3ra  +  *  2U1/  Ifu^tt  . 

2 

.low  the  rlsht-hand  member  converges  strongly  to  0  in  L  (p; 


since  the  Incident  field  u*  is  fixed,  while  the  compactness 
of  the  operators  D  and  3  insure  that  2  D 4  and 

ill 

3  «r  ^  -3  a  .  jut  then,  the  relation  (12;  can  be 

m  •  m  i 

rewritten  as 

(13)  *\)  =  2U1/  11  u  41  -  3  or  if>  -  if 

m  in  m  ,a  m 

which  shows  that  .  Moreover,  this  strong  convergence 

of  the  ^  is  enough  to  insure  that  the  products 
in  L2(0)  and,  consequently ,  that  the  function  ^satisfies 

(14)  (I  +  5  cr  +  2)^  =  0. 

Mow,  consider  the  second  boundary  integral  equation 

(15)  (-  cri  +  K*  a  +  D  )  u  =  2  fcuVin 

n 

ana  look  at  the  sequence 

(16)  D  ^  *  2(  **  uVit)  tlu  W  1  -  (  -  O'  +  n*  a  )  _  . 

n  m  m  s  a  1  a 

3ince  the  *1*  converge  to  ^  strongly  and  the  oroducts  C r  it 
m  *  a  1  m 

converge  to  weakly  in  L2(f),  the  functions  Dn^  defined 
by  (16)  must  converge  weakly  in  this  space  to  . 

This  fact,  together  with  the  estimates  which  appear  in  i.  3  ]» 
shows  that  the  function  f  ,  as  a  solution  of  the  first  equation 
(12),  must  lie  in  the  domain  of  the  operator  Dn  and  must  satisfy 
the  homogeneous  equation 

(17)  (.»I  +  l*w  +  Dn}^  =»  0. 

ilence  the  function  ,  as  the  strong  limit  of  elements  of  norm 
1,  must  be  a  non-trivial  solution  to  the  homogeneous  system  of 
boundary  integral  equations,  which  contradicts  the  unique 
solvability  of  the  system  (9). 

AC:: ;0.i  L23  O.c.  MT :  This  work  was  supported  by  the  U.5.  Air  Force 
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